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We present a quantitative semi-classical theory for the non-equilibrium dynamics of transverse 
Ising chains after quantum quenches, in particular sudden changes of the transverse field strength. 
We obtain accurate predictions for the quench dependent relaxation times and correlation lengths, 
and also about the recurrence times and quasi-periodicity of time dependent correlations in finite 
systems with open or periodic boundary conditions. We compare the quantitative predictions of our 
semi-classical theory (local magnetization, equal time bulk-bulk and surface-to-bulk correlations, 
and bulk autocorrelations) with the results from exact free fermion calculations and discuss the 
range of applicability of the semi-classical theory and possible generalizations and extensions. 



I. INTRODUCTION 

The non-equilibrium quantum relaxation in many- 
body systems has gained increased interest over the re- 
cent years, not least because trapped cold atom systems 
made its experimental study possible. In principle one 
asks for the fate of an initial state that is not an eigen- 
state of the Hamiltonian under the time evolution ac- 
cording to the Schrodinger equation. A straightforward 
method to prepare such an initial state is the instanta- 
neous change of a global or local parameter of the system 
like an external field or the interaction strength, denoted 
as a quantum quench or simply quench. Important issues 
of interest are then: 1) Is there an asymptotic stationary 
state, what are its characteristics, is it describable by a 
general Gibbs ensemble (i.e. does the system thermalize 
after a quench)? 2) What are the characteristics of the 
dynamical evolution of order, correlations and quantum 
entanglement in the system? 

The first theoretical studies of quenches in quantum 
many body systems were performed for the quantum 
XY and quantum Ising spin chains^"— . Spectacular ex- 
perimental results^ triggered an intensive research on 
quantum quenches in various systems such as ID Bose 
systems^, the quantum sine-Gordon model^, Luttinger 
liquids^ and othersiS,. Besides studies on specific models 
there are also ficld-thcorctical investigations, in which re- 
lation with boundary critical phenomena and conformal 
field-theory are utilizedii"— . Progress in understanding 
thermalization, or absence thereof, in a particularly well 
studied integrable model, the transverse Ising chain, has 
been achieved iniiii^. The concept of an effective tem- 
perature depending on the quench parameters is useful 
to parametrize the relaxation time and correlation length 
determining the spin correlations after a global quench. 
But actually each excitation mode has its own thermal- 
ization temperaturoi^, implying that the system never 
thermalizes after a quench. 

For the transverse Ising chain in thermal equilibrium 
Sachdev and Youngii introduced a semi-classical de- 



scription of the equilibrium quantum relaxation in terms 
of ballistically moving quasi-particles. This description 
turned out to be surprisingly accurate in predicting the 
temperature dependence of relaxation time, correlation 
length and scaling forms in the ferromagnetic and para- 
magnetic phase. 

For global quantum quenches a picture of ballistically 
moving quasi-particles spontaneously created after the 
quench has been usediiii^ to explain several features of 
the time-evolution of different quantities, in particular 
that of the entanglement entropyi2i2£. This picture has 
also been used to interpret results of exact calculations 
obtained with the free fermion techniqu o^^'^^ or field the- 
ory (at the critical point) ^^i^^ . 

Obviously it would be desirable to have a quantita- 
tive semi-classical theory for the non-equilibrium dy- 
namics after quantum quenches, too. This is what we 
will present in this paper for global quenches, for local 
quenches a brief account has been given by us recently 
in2i. Here we present the quantitative analog of the semi- 
classical theory for equilibrium quantum relaxation of 
transverse Ising chainaii and generalize it to the non- 
equilibrium dynamics in finite systems. By this we will 
not only obtain accurate predictions for the relaxation 
times and correlation lengths, but also about the recur- 
rence times and quasi-periodicity of time dependent cor- 
relations in finite systems with open or periodic bound- 
ary conditions. Since in experimental set-ups of quantum 
quenches, as for instance cold atom systems, the number 
of particles is rather restricted and far away from the in- 
finite system size limit, the understanding of finite size 
effects in non-equilibrium quantum relaxation is impor- 
tant and may be, as we will show often be drastic. 

The paper is organized as follows: After the model def- 
inition in the next section we present the semi-classical 
theory for the non-equilibrium dynamics of the trans- 
verse Ising chain after a quench. Then we derive the 
semi-classical formula for the local magnetization, equal 
time bulk-bulk and surface-to-bulk correlations, and bulk 
autocorrelations and compare the predictions with the 
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results from exact free fcrmion calculations. Finally we 
discuss the range of applicability of the semi-classical the- 
ory and possible generalizations and extensions. 

II. MODEL 

The system we consider in this paper is the quantum 
Ising chain defined by the Hamiltonian^^: 

1=1 1=1 

in terms of the Pauli- matrices af'^ at site I. In ^ the 
chain has a finite length L and open boundaries, later 
we will also discuss periodic boundary conditions. We 
consider global quenches in which the transverse field 
strength is suddenly changed from ho ioi t < to h ^ ho 
for t > 0. For < < the system is in equilibrium, which 
means it is in its ground state j'Z'o) and which we denote 
as its initial state. After the quench, for t > 0, the state 
jtf'o) evolves according to the new Hamiltonian: 

\%it)) =expi~im)\%) . (2) 

Similarly we have for the time-evolution of an operator: 
ai{t) = e'xp{—it'H)ai expiitH). 

We consider the general, time and space dependent 
correlation function: 

C{r,M;r2M) = {M<,(ti)<,(t2)W^) , (3) 

and study its behavior in special circumstances. The 
autocorrelation function is obtained for ri = r2 = r, 
which is denoted as Gr(ti,i2), whereas for ti = t2 ^ t 
we have the equal-time correlation function. This latter 
quantity for large separation behaves as: C(ri, <; r2,t) = 
Ct{ri,r2) = mr^{t)mr^{t), where mr{t) is the local mag- 
netization. In the initial state (and in the thermody- 
namic limit, L — > oo) for ho < he = 1 there is a finite 
magnetization, niriO) > 0, whereas for ho > 1 one has 
mr{0) ~ 0{1/L). 

A. Free fermion representation 

The Hamiltonian in Eq.(IT]) can be expressed in terms 
of free fermion creation, 7]^, and annihilation operators, 

H^Y.^H{p){vlVp~ll'^) , (4) 
p 

where the energy of modes is given by 

£h{p) ^ \/{h~ cosp)2 + sin'^p . (5) 

The quasi-momenta, p, has L quasi-equidistant values 
in the interval: < p < tt for free boundary condi- 
tions, whereas for closed chains these are restricted to 



\p\ < TT. Time-evolution of the fermion operators are 
'7p(0 = and f?fc(t) = e~'*^'=77fc from which one can 

obtain the time-evolution of the spin operators. The cor- 
relation functions in the fermion representations are ex- 
pressed in terms of Pfaffians, which arc then calculated 
as the square-root of the determinant of the correspond- 
ing antisymmetric matrix, which has the elements of the 
Pfaffian above the diagonal. For free boundary condi- 
tions these determinants have a dimension: 2(ri + r2). 
Following Yang2i, the local magnetizations can be cal- 
culated in the form of an off-diagonal matrix-element: 
mr{t) = {%\cT^ {t)\^i) , where denotes the first ex- 
cited state for t < 0. Its numerical calculation necessities 
the solution of a 2r x 2r determinant. 



III. SEMI-CLASSICAL THEORY 

In the absence of the transverse field in Eq.((T]), h ~ 0, 
the system is, in the transformed basis af ^ erf identical 
with the classical Ising spin chain. The ground state is 
two-fold degenerate and given by: |!Z'o) = I + + + ••• +) 

and I'Fq) = 1 ■ ■ ■ ^) and the first excited states are 

the (L — l)-fold degenerate given by the single kink states 

\n) = \ + + ■ ■ ■ + -\ • • ), where n denotes the kink 

position. Switching on a small transverse field, h > 0, the 
low-lying excitations arc, in first order order degenerate 
perturbation theory, superpositions of these single-kink 
states X^n'^"!'^) with excitation energy eh{p)- The ac- 
tual perturbation calculation yields a,i = ■\/2/i sin(pn) , 
with £h{p) = 1 — ft. cosp, where p has i — 1 discrete values 
in the same region as given below Eq.(l5|). Thus the low 
lying excitations of H are Fourier transforms of localized 
single kink states, similar to the eigenstates of the Hamil- 
tonian for free particles in a box of length L. Analogously 
freely moving single kinks are therefore wave packets of 
the aforementioned low lying excitations. Their energy 
agrees to to leading order in h with the free-fermion en- 
ergies in Eq.([5]) and they move ballistically with constant 
velocity ±Vp given by 

Vp = ^ = . 6) 

dp ep 

Ballistically moving kinks are then the (fermionic) 
quasi-particles (QPs) which we use in the following to 
formulate a semi-classical theory of the quantum quench 
dynamics of the transverse Ising model. Since by defi- 
nition these QPs are well-defined at small fields in the 
ferromagnetic phase, the theory is expected to be appli- 
cable for quenches in the ferromagnetic phase. It will 
turn out that it actually holds in the whole ferromag- 
netic region not too close to the critical point {h = 1). In 
the paramagnetic phase one can start with the /i — > oo 
ground state to introduce an analogous QP-concept in- 
volving individual spin flips instead of kinksii but the 
same dispersion relation and velocity We will 
mention the necessary modifications below. 
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Immediately after the quench the time-dependent state 
of the system in Eq.(I2]), which for small h and for small 
t is given by 

\Mt)) - cxp (^-tthJ2^f^ 1*^0) 

= Yl [cos{th) + ism{th)af] \%) . (7) 

(where we use the convention, erf o af , as before) . This 
indicates that by the action of the af operators initially 
single spins are flipped and thus pairs of kinks are created 
at each lattice point, which then move ballistically with 
a speed v ^ h. The maximum velocity is Vmax ~ h for 
small h. 

In a translationally invariant system the creation prob- 
ability of QP-s is uniform and will be denoted by 
fp{h{),h). For open boundary conditions there will be 
corrections to a uniform creation probability close the 
boundaries, which are negligible for sufficiently large sys- 
tem sizes. In an equilibrated system that is thermalized 
at temperature T, this would be f°^{ha,h) = e"'^"/-^. 
For zero temperature quantum relaxation fp{ho, h) is the 
probability with which the modes with momentum num- 
ber p are occupied in the initial state |5'o), i-e. 

(8) 



fpiho,h) = (*ohp?7p|*o)- 



In a finite system with open boundaries a QP with 
momentum p moves uniformly with velocity Vp until it 
reaches one of the boundaries, where it is reflected and 
moves with velocity ~Vp thereafter, and so forth. The 
trajectory of the kink is periodic in time, after a time 



2Tp with 



Tp = L/vp 



(9) 



(including a reflection at the right and left boundary) it 
returns to the starting point xq with the initial direc- 
tion and velocity Vp, see Fig. [1] Due to conservation of 
momenta after a global quench QPs emerge pairwise at 
random positions with velocities +Vp and —Vp, as indi- 
cated in Fig. [T]for three QP pairs. 

For a given QP pair created (at i = 0) at position 
xo e [0, L] let xi{t) be the position of the initially right- 
moving QP (i.e with initial velocity Vp) at time t and 
X2(t) be the position of the initially left-moving one (i.e. 
with initial velocity —Vp)- Deflne ta as the time when the 
left-moving particle reaches the left wall the first time and 
tb as the time when the right-moving particle reaches the 
right wall the first time 



tb 



xo/vp 

{L - xo)/vp 



(10) 



Then for t < T„ 
xi{t) = 

X2{t) = 



Xq + Vpt 
2L — Xq — Vpt 



Xq - 
Vpt 



Vpt 
- Xq 



for t < tb 
for tb <t<Tp 

for t < ta 

for ta <t <Tp 




FIG. 1: (Color online) Typical semi-classical contribution to 
the correlation function C{ri,ti;r2,t2)- Note the 6 trajec- 
tories of the 3 QP-pairs intersect the line (ri,ti;r2,t2) five, 
i.e. an odd number, of times, which implies that a^^{ti) and 
0-^2(^2) have opposite orientation. Equivalently one can say 
that the trajectories of the red and the green QP pair inter- 
sect (ri, ti; r2,t2) an even number of times (and thus do not 
contribute) and the trajectory of the blue QP pair an odd 
number of times. 



At t = Tp the two QPs meet at x — L — xq. For Tp < 
t < 2Tp the trajectories are defined accordingly (see Fig. 
[2]), and for t > 2Tp one notes that xi and X2 are 2Tp- 
periodic. 

Since QPs represent kinks or domain walls, changes 
sign each time a QP passes. Therefore the correlation 
function in Eq.Q can be evaluated in terms of classical 
particles moving according to by using a similar rea- 
soning as in equilibriuroii, the difference being that here 
1) QP trajectories can intersect the line (ri, <i; r2, ^2) sev- 
eral times, 2) QP trajectories come always in pairs with 
a common off-spring at i = 0, and 3) the occupation 
number of QPs is not thermal. 

If a QP trajectory intersects the line (ri, ii; r2, ^2) an 
odd number of times, the spins at (ri,ti) and (r2,t2) 
have the opposite orientations (i.e. crfi(^i) = — crr2(*2)), 
which contributes to the decay of the correlation between 
<Jr^{ti) and a^^{t2), see Fig. [TJ If the two trajectories 
pass an even number of times, the spins have the same 
orientation, as if the trajectories did not pass the line 
iri,ti;r2,t2) at all. Let Q{ri,ti;r2,t2) be the probabil- 
ity, that the QPs, which have started from the same site, 
have passed the line {ri,ti; r2,t2) a total odd number of 
times. Then the probability that for a given set of n sites 
the kinks have passed (for each site total odd times) this 
line is: Q"(l — Q)^^". Summing over all possibilities we 
have: 



C{ri,ti]r2,t2) 

Ceq(ri,r2) 



|](-l)"g"(l-Q)^-" 



LI 



n\{L~n)\ 



(1 - 20)^ « e-2Q('-i.*i^'-2,t2)L ^ ^2) 



(11) 



where Coq(ri,r2) is the equilibrium correlation function 
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in the initial state and in the last step we have used 
that the probability, Q{ri,ti;r2, ^2) is small. To calculate 
Q one should average over the QPs with momenta p G 
[— tTjTt], or equivalently one can average over QP-pairs 
which is restricted to p € [0, tt]. In this second method 
we have the expression: 

1 r 

Q{ri,ti;r2,t2) = — j dp .fp{ho,h) ■ qp{ri,ti;r2,t2) 

(13) 

in terms of the occupation probability (see Eq.®) and 
the passing probability, qp{ri,ti;r2,t2)- This latter 
quantity measures the probability that the two trajec- 
tories xi(t) and X2{t) of any QP-pair with momentum 
p intersect the line ("Ti, ii; ^2, ^2) together an odd num- 
ber of times. The same probability for a given QP- 
pair which is emitted at site xq € [0,L] is denoted by 
qp{xo\ri,ti;r2,t2)- If we assume that the generation of 
QPs at the quench is homogeneous in space then we ob- 
tain: 

qp{ri,ti;r2,t2) ^ Y [ dxo qp{xo\ri,ti;r2,t2) ■ (14) 
^ Jo 

In most cases of interest (see below) it is possi- 
ble to provide an analytical form for the function 
qp(j'i,ti;r2,t2)- If not, the number of intersections 
can straightforwardly be determined numerically and 
averaged over xq, yielding (3'p(f 1, ii; ^2, ^2) and thus 
Q{ri,ti;r2,t2) in Eq. p^ and the correlation function 
in Eq.®. 



t<l/v 




FIG. 2: Left: Typical semi-classical contribution to the time 
dependence of the local magnetization mi{t). Full lines are 
quasi-particles or kinks moving with velocity Vp through the 
chain. The ± signs denote the sign of the spin at site I. 
Right: Sketch of the trajectories of kink pairs that flip the 
spin at position I exactly once for times t < Tp/2. Kink pairs 
with initial position xq outside the marked region either do 
not flip the spin at / (since they do not reach the position I 
within time t) or they flip it twice, qp is the fraction of the 
marked intervals on the t — 0-axis. 



IV. LOCAL MAGNETIZATION 

The time-dependent local magnetization at a site I 
(here we consider I < L/2) can be formally expressed 
to a correlation between a spin that is fixed at time t = 
(to, say, af = +1) and the same spin at later times t, i.e. 
mi{t) = m^'i • C\^.(^t=o)=+{l,0;l,t). Then, with Eq.® 

m,(t) =m;='i-e-29(*.')^ (15) 

with q{t,l) = Q|o-^(t=o)=+i(^ 0; /, t), which is with 
Eq.dni) 

<litJ)^^ f dp fp{ho,h)qp{t,l) (16) 

where 

9p(*'0==^y" dxoqpixo,tJ) (17) 

as in Ea. ([T^ . To calculate qp one concentrates first on 
times t < Tp/2 = L/2vp. Now qp is just the fraction of 
possible initial positions from which kink pairs can start 
with velocity +Vp and —Vp that flip the spin at position I 
exactly once. This region is marked in the sketch of Fig. 



m One sees that for t < l/vp one gets qp = 2vt/L, and 
for l/vp < t < Tp/2 one gets qp = 21/L, independent of 
time. 

For Tp/2 < t < Tp one observes that a kink pair that 
started (at i = 0) at position xq re-unites after a time 
t = Tp at position L — xq. Since the origins of kink pairs 
are distributed uniformly over the chain the probability 
qp{t,l) is Tp-periodic (n.b.: the kink trajectories them- 
selves are only 2rp-periodic). Moreover qp{t,l) is sym- 
metric with respect to time inversion since it is symmet- 
ric under the QP velocity inversion, qp{-~t,l) ~ qp{t,l), 
therefore qp{Tp ~ t,l) = qp{t,l). Defining the reflection 
times ti = l/vp and t2 = Tp — ti one then has for the 
period < i < Tp for ^ < L/2 

{2vpt/L for t < ti 

21/L for h<t<t2 (18) 

2 - 2vpt/L for t2 <t <Tp 

For / > L/2 one uses the symmetry qp{t, I) = qp{t, L — I) 
and for t > Tp one makes use of the Tp-periodicity of 
qp{t,l): 

qp{t + nTp,l)=qp{t,l), (n=l,2,...). (19) 

Although qp(t,l) is Tp-periodic, q(t,l) is not periodic, 
since all QPs have different speed. Nevertheless the max- 
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FIG. 3: Semi-classical prediction for the local magnetization 
mi{t) quench. Here L = 1024, ho ^ 0, h = 0.2. The (quasi)- 
periodicity ^0} is Tperiod = L/h = 5120. 



imum speed fmax = h + 0{h?) determines the onset 
of magnetization reconstruction and therefore a quasi- 
periodicity of q{t,l) and concomitantly mi{t), whose 
(quasi)-period is then expected to be 

^period = i/Wmax ~ L/h (20) 

With Eqs.dini), (HSI) and dH]) one obtains mi{t) via nu- 
merical integration (or summation over the discrete p- 
values for a lattice of finite size L). 

For an actual calculation one needs to know the occu- 
pation probability fp{hQ, h) in Eq.Q, which can be cal- 
culated numerically in a straightforward manner using 
the free fermion technique. Since for large system sizes 
the occupation probability is not expected to depend 
strongly on the boundary condition we will use in the 
following the expression for fp{ho, h) for periodic bound- 
ary conditions, which can be given in analytical form as 
shown in Appendix A, see Eg. ([55]) . In Fig. [3] the pre- 
diction of the semi-classical computation is shown. One 
observes the predicted quasi-periodicity for finite lattices 
and the expected exponential decays in / and t as dis- 
cussed below. 

In Fig.|3]we compare the semi-classical prediction with 
the exact results obtained with the free fermion technique 
and find that the agreement is remarkably good. We ob- 
serve that small deviations occur in the bulk {I ^ L/2) 
for t > Tporiod/2, which is when the first QP reflec- 
tions are involved in the dynamical evolution of mi{t). 
For sites close to the boundary, i.e. small I we observe 
small deviations already in the plateau region (c.f. also 
the surface-to-bulk correlations discussed further below) . 
Here a spatially inhomogeneous QP creation probabilities 
would have the most significant effect, which is negligible 
for bulk spins. 

When the system after the quench would be ther- 
malized at some effective temperature Tce{hQ, h)^ this 



FIG. 4: Relaxation of the local magnetization, logm;(f), at 
different positions in a L = 256 chain with free ends after 
a quench with parameters ho = 0.0, h = 0.2 and L = 256. 
A Exact (free fermion calculation). B Semi-classical predic- 
tion (|15() with the passing probability (|18p and the occupation 
probability (|55|) . C Comparison between exact and QP cal- 
culation for mi{t) for L — 256, / = 128 for a quench from 
fto = to ^ = 0.1. D Semi-classical prediction using a ther- 
mal occupation number probability in Eq. (|21P with an effec- 
tive temperature, Tcit, see the text. 

would imply that the occupation probability is 

fp{ho, h) = e-^''('')/^'--«(''°''') . (21) 

The effective temperature is determined from the con- 
dition, that the relaxation time in equilibrium tt(/i,T) 
(with transverse field h and temperature T) is the same 
as in quantum relaxation at T = but after a quench 
from ho to h. In the limit T < A(/i), A(/i) being the 

gap of the system we haveii: tt(/i, T) sa '^^'^^^ ^ which 
should be compared with t(/io, ft.), which for small h and 
/iQ is given in Eg . ([55)) . The result of the semi-classical 
calculation using (|2ip is also shown in Fig. Hp and com- 
pared with the exact data for a quench form /iq = to 
h = 0.2. One sees that, using the proper effective temper- 
ature the initial exponential decay agrees perfectly, but 
soon as the first reflections are involved, large deviations 
occur. An effective temperature can describe the initial 
relaxation well because essentially it is a fit parameter 
for the initial exponential decay, which stops after some 
(Z-dependent) time in a finite system. 

In the infinite system size limit L — > oo the time ^2 
in (|18p is infinite for all momenta p. Thus the (quasi)- 
periodicity of mi{t) is lost and the functional form of 
771/ (t) as predicted by (HH]), (US]), and (HH) is 

mi{t) = exp^-t-^y dpvp fp{ho,h) 9{l - Vpt)^ 

■ exp(^-l-^J^ dp fp{ho,h)0{vpt- 1)^22) 
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FIG. 5: Relaxation of the local magnetization, m;(t) after 
quenches into the disordered phase - from the ordered phase 
{ho = 0.5, h = 1.5) (A) exact, (B) semi-classical, and from 
the disordered phase (ha — 1.5, h = 2.0) (C) exact, (D) 
semi-classical. The legend of A and B holds also for C and D, 
the system size is L = 256. 



which defines, in analogy to^ the quench specific length 
and time scales 



mag 



lJhQ,h) = 



dpvp fp{ho,h) 



Cmag(^0,/l) = - / dpfp{ho,h) 



(23) 



In the small h and ho limit these are calculated in 
Eqs.([39l) and respectively. 

In a finite system there is a quasi-periodicity and the 
magnetization after the first relaxation period is recon- 
structed. Due to the p-dcpendence of the velocity of the 
QPs in the reconstruction regime the rate of exponential 
increase of the magnetization, Tj^^^g, is increasing in time. 
Its maximal value is reached ait = Tperiod, which is given 
by: 



1 2 


'f ' 








Jo 



dpVp fp{ho,h) 



h{h — ho 



,9^3- 



127r 



(24) 



where the second expression is valid in the small h and 
limit. One can see that r(/io,/i) < T'{hQ,h), thus 
the reconstruction is slower than the relaxation. While 
QPs with large energy and high velocity contribute to 
the reconstruction, the other QPs with smaller energy 
and lower velocity still reduce the magnetization. These 
processes with opposite effect are responsible for the de- 
cay of the amplitude of the quasi-periodic oscillations of 
the profile, see FigsIS] andH) 

After quenches into the disordered phase {h > 1) the 
relaxation (and recurrent) dynamics of the longitudi- 
nal magnetization is superposed by oscillations from the 



ground state correlationaii and one has to replace rn^'^ 
in ^ by 



K{tA) 



(25) 



where K{x) is the modified Bessel function. The results 
for the corresponding QP calculation and comparison 
with the exact data are shown in Fig. [5l One observes 
again that the relaxation and recurrent dynamics is well 
described by the semi-classical picture also for quenches 
into the paramagnetic (disordered) phase. The super- 
posed oscillations have a slightly larger amplitude and 
frequency. Note also that for quenches from the para- 
magnetic phase (Fig. [S|) C, D) the equilibrium profiles 
771°'^) shifts the curves for mi{t) downwards for increasing 
I, since in the paramagnetic phase the surface magnetiza- 
tion is larger than the bulk magnetization in a finite chain 
(both vanishing only in the infinite system size limit). 



V. CORRELATION FUNCTIONS 

As mentioned before it is possible to perform the 
semi-classical calculation for the two-spin correlations 
C(ri, ti;r2, ^2) for any pair of sites ri, r2 and any pair of 
times ti, t2 with the formulas ([TB]). and (|14p . Here 
we want to focus on the time dependence of equal time 
correlations between spins separated by a distance r and 
arranged symmetrically within the bulk, i.e. we consider 

Ctir) = C{L/2 - r/2, t; L/2 + r/2, t) , (26) 

for quenches within the ordered phase (/i < 1), which is, 
within the semi-classical theory given by: 

Ct{r) = Coq(r) • exp dp fp{ho, h) ■ qp{t, r) 

(27) 

As sketched in Fig. |6]the function qp{r,t) for Ct(r) is 
Tp/2-periodic and for the period < ^ < Tp/2 given by 
for r < L/2 



4vpt/L 
2r/L 

2 - Aupt/L 



for 
for 
for 



t < ti 

h<t<t2 

t2<t < Tp/2 



(28) 



with ti = r/2vp and t2 — Tp/2 — ti. (For r > L/2 one 
should replace in the above formulas r to L — r.) Note 
that the relevant times occurring in this expression are all 
multiplied with a factor 1 /2 as compared to those deter- 
mining Qp for the local magnetization (jl8p . In particular 
qp(t,r) is rp/2-periodic (in contrast to the Tp periodic- 
ity of Qp for the local magnetization): qp(t + nTp/2, r) = 
qp{t, r) for n = 1, 2, 3, . . .. As a result the (quasi)-period 
of Ct{r) for fixed r is one half of the (quasi)-period of the 
local magnetization mi{t) 



'^period — L/2va 



L/2h 



(29) 
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(A) 



t < r/2v t 



q = 4vt / L 2vt 



t > r/2v t 
t < Tp/4 



FIG. 6: (Color onlne) Semi-classical contributions to the 
equal time correlation function Ct(r) = C(L/2 — r/2, t; L/2 + 
r/2,t). Left: Sketch of the trajectories of kink pairs that 
reverse the orientation of the spins at ri and r2 for times t < 
Tp/4:. Kink pairs with initial position xo outside the marked 
region either do not intersect the line (ri = L/2 — r/2, t; r2 ~ 
L/2 + r/2, t) (red) (since they do not reach the red line within 
the time t) or they flip it twice. Qp is the fraction of the marked 
intervals on the t — 0-axis. Right: Sketch of the additional 
symmetry of qp{t) that reduce its periodicity from Tp to Tp/2: 
For each QP pair created at position xq intersection the red 
line at time t < Tp/A there is a QP pair created at position 
L — Xo that intersects the red line at time Tp/2 — t. Hence 
Ipit) = <lp{Tp/2 - t) for t < Tp/2. At t = Tp/2 the QP pair 
created at xo meets again at L — xq and the one created at 
L — Xo meets again at xo, which implies after averaging over 
initial position that qp{t + Tp/2) = qp{t). 



With (p8)) and fp{ho,h) from appendix A the semi- 
classical calculation can be performed, results and the 
comparison with exact data are shown in Fig. [7] and [51 

In[n]wc show the semi-classical prediction for Ct(r) for 
larger system sizes and long times, scaled by the (quasi)- 
period T^^^-^^^, ((^ . which demonstrates the persistence 
of the recurrence for very long times in finite systems. 
Note that the recurrence amplitude decreases with in- 
creasing system size and vanishes completely for L — >■ oo. 

In the infinite system size limit L — >■ oo the time t2 
in (jlSp is infinite for all monrenta p. Thus the (quasi)- 
periodicity of C't{r) is lost and the functional form of 
Ct(r) as predicted by ^ and is 



Ct{r) = Ccq(r) exp dpvpfp 0{l - Vpt) 



exp [ -r • - y dp fp e{vpt - I) ) (30) 



with fp = fp{hQ,h). This agrees to first order in fp to 
the prediction o&, where fp is replaced by —1/2 log(l — 
2/p) = fj^ + 0[fl) (see appendix A). Eq. ^ defines the 
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FIG. 7: Equal time correlation function Ct(r) for fixed r as a 
function of time t after the quench: Comparison between the 
exact result (left) and the semi-classical prediction (right). 
L — 256, ho = 0., h — 1.5 (A) exact, (B) semi-classical; 
ho = 0.3, h = 0.5 (C) exact, (D) semi-classical; The legend 
of A and B holds also for C and D. 



1 



Q g _ 'l| ^lLllllJlJJJJj:j.LULlLnJlJjJj;jJ.LLLLUllJlJJJJJJJ.LLLLUllJJJJJJ. i .L i .L i .UUjlJlJJJJ. i . i . i .LLLUlllJ_, - 



0.8 

_ 0.7 
^ 0.6 

0.5 



0.4 



t=10 
t=20 
t=40 
t=60 



20 40 60 80 100 
r 



FIG. 8: Equal time correlation function Ct(r) for fixed time 
t after the quench as a function of distance r Comparison be- 
tween the exact result (points) and the QP calculation (black 
broken lines). L = 256, ho ~ 0, and h = 0.25. 



quench specific length and tinre scales 



dpvp fpiho^h) 



Sc^(^o,/i) - / dpfp{ho,h) 



(31) 



Note that Tc = Tmag/S and = 6nag. For a small ho and 
h this yields to leading order 

C '{ho,h)^ — / dfcsm k= (32) 
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FIG. 9: QP prediction for the equal time correlation function 
Ct{r) plotted against time after the quench scaled with the 
(quasi)-period Tj^i-iod = L/2h for different fields h {ho = 0). 
The three curves for each field h correspond to different sys- 
tem sizes: L — 256, 512, and f 024 (from the top curve to the 
bottom curve). 



FIG. fO: Surface-to-bulk correlation function Cf"''(r) for 
fixed time r as a function of the time after the quench. Com- 
parison between the exact result (points) and the QP calcula- 
tion (black broken lines). L = 256, ho = 0.75, and h — 0.25. 



B. Autocorrelations 



Surface-to-bulk correlation 



The surfacc-to-bulk correlation function C™''^(r) 
C(0, t; r, t) is within semi-classical theory given by 



Cr'(r)=C,7f(r).exp(-- / dp ^ 



(33) 

with fp = fp{hQ,h). Similar considerations that lead to 
([T8|) and (|28p yield an analytical expression for (/^"'^^(t, r), 
which is equivalent to qp{t, I) for the local magnetization 
in (1181). however with I = r: 



qp{t,l = r) 



Consequently 



mi=r{t) 

CO 



(34) 



(35) 



which implies that the surface-bulk correlation is domi- 
nated by the relaxation of the magnetization at the bulk 
site and that it is Tp-periodic (in contrast to Cf (r), which 
is Tp/2-pcriodic. 

In Fig. [To] we show a comparison of this semi-classical 
result with the exact data. We observe that small devi- 
ations occur in the plateau region for small distances r 
where the bulk-correlations still agree very well with the 
semi-classical prediction. Since for small r both sites in 
the surface-to-bulk correlation function are close to the 
boundary a spatially inhomogeneous QP creation proba- 
bilities would have the most significant effects here. 



The autocorrelation function 

Gi{t) = C{l,Q;l,t) 



(36) 



is (up to an extra factor rnp-) identical to the time- 
dependent local magnetization mi{t): 



Gl{t) — Geq gq^ 



(37) 



For I = L/2 (bulk autocorrelation) in the limit L — >■ oo 
the QP prediction is 



Gl/2(*) oc exp 



dp fp{ho,h) ■ Vpt 



-t/Tauto 



with the relaxation time 



' mag J 



(38) 

ea. (|23)) . which cor- 



responds to the leading order of the result from Calabrese 
et aJi^(see appendix A), in which fp is again replaced by 
-1/2 log(l - 2/p) = fp + 0{fl). For a small ho and h 
([23)) yields to leading order 



27r 



^_,^ h{h hof f\ksm'k^hih~ho)''— . (39) 
Jo 37r 



This has already been found numerically i: 



as 
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VI. PERIODIC BOUNDARY CONDITIONS 

In a chain with periodic boundary conditions instead of 
the open boundaries that we considered so far, one has 
to replace the QP trajectories ([TT|) by the appropriate 
expressions: 

Xi{t) = (.To + Vpt) mod L 

X2{t) = {xq — Vpt) mod L (40) 

where the modulo operation is defined in the obvious 
manner: Shift the real number Xi by multiples of L such 
that it lays in the interval [0, L]. With this the evaluation 
of the local magnetization and correlation functions is 
straight-forward. 

The chain with periodic boundary conditions is trans- 
lationally invariant, therefore the equal time correla- 
tion Cf'^''^'{r) — C^-^'^ {ri,ti;ri + r, ^2) is indepen- 
dent of ri. One sees immediately that the expression 
for gP '^ '= (t,r) is identical to q^{t,r) in Eq. ([Ml), and 

therefore Cf' '^{r) is identical, up to prefactors from 
the ground state or equilibrium correlation function, to 
C°P™(L/2 - r/2, t; L/2 + r/2, t): 

fjp.h.c.f \ 

C?''-'ir) = -^§^-Cr''{r) (41) 

Ooq {r) 

It should be noted that this relation only holds for the 
symmetric correlation function C°^°'^{r) = C°p°"(L/2 — 
r/2,t; L/2 + r/2, t). 

The local magnetization mi{i) is independent of the 
site / in a system with periodic boundary conditions, as 
is the QP passing probability gP '' '= (t, = (?P '' '= (t). We 
find 

n.b.c./,N _ / 2vpt/L for t < Tp/2 

" \ 2 - 2vpt/L for Tp/2 <t<Tp ^ > 

For t > Tp one uses the Tp-periodicity q^-^ '^-{t + nTp) = 
qP-^ '^ it), (n = 1,2, . . .). With eq. ^ and (HH) the local 
magnetization is then given by 

p.h.c.(-j-\ f 2 r 

pb, = exp — / dpfp Vpt e{sm{2nt/Tp)) 
rrieq ■ ■ V Jo 

+ dp fp{L-Vpt)d{-sin{2TTt/Tp))^ (43) 

In the infinite system size limit this yields 

mP ''-'= (t) oc exp (^-t^ dp fp Vp^ cx 6-*/^"""= , (44) 

which agrees to first order in fp with the prediction oi^. 
As for open boundary conditions the autoscorrelation 
function G^^^■'^■{t) is given, up to prefactors, by the same 
expression as the local magnetization (j43|) . 



VII. DISCUSSION 

We have formulated a semi-classical theory for the non- 
equilibrium quantum relaxation of the transverse Ising 
chain after a global quench via an instantaneous change 
of the transverse field. It is applicable to systems of finite 
and infinite length and describes properly the relaxation 
dynamics as well as the recurrence / reconstruction prop- 
erties of dynamical correlations in finite systems. For in- 
finite systems our theory agrees to lowest order with a re- 
cent prediction by Calabrese et ali^. Our results indicate 
that the global quantum quench induces a unique length- 
scale, ^, and a unique time-scale, t, in the system, both 
dependent upon the quench parameters, ho and h. These 
characteristic scales appear also in half-infinite geometry 
and in finite systems, provided the length of the system 
is larger than ^. In a finite system this semi-classical the- 
ory not only explains the recurrence and reconstruction 
properties of the local magnetization^^, but describes the 
dynamical behavior quantitatively. 

The semi-classical theory can be used to define an ef- 
fective temperature for the quantum relaxation process. 
If we compare the expressions obtained by Sachdev and 
Youngii for the correlation length and the relaxation 
time in equilibrium at finite temperatures with our re- 
sults for zero temperature quantum quenches one obtains 
a node-dependent effective temperature, Tcff(p), defined 
by the condition: 

/p(/io,/i) = exp(^-g^) . (45) 

This relation agrees to first order in fp with the prediction 
oSi^ (i.e. for small effective temperatures or small differ- 
ences \h — hoi). Ini^ the Boltzmann-factor on the r.h.s. 
of Ea. (|15)) is replaced by the Fermi- function with zero 
chemical potential, as shown in the Appendix in Ea. (|58p . 
thus replacing classical kinks simply by free fermions. 

It is interesting to notice an analogous expression for 
the time-evolution of the entanglement entropy, S{t), 
measured after the quench between two semi-infinite 
parts of the system, say A and B. The analytical re- 
sult by Fagotti and Calabrese^ can be written into the 
form: 

1 r 

S{t)=t- dpVpSp{h,ho) , (46) 
Ti" Jo 

with 

3p{h,ho) = -(1 - /p)ln(l - fp) - fplnfp , (47) 

being the entropy of the fermionic mode with occupation 
number fp{ho,h). In the semi-classical theory this ex- 
pression can be interpreted as the result of ballistically 
moving QP pairs, which are created say at ^ at t = 
and one of them is reaching B before the actual time, t. 
Each of these QPs brings an entropic contribution as a 
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free fermion. It would be interesting to see if the rela- 
tions in Eas. ([55)) and pS)) are valid for another integrable 
quantum spin systems, too. 

The semi-classical approach is accurate, if the occu- 
pation probability, fp{ho,h), is small, which is valid if 
the initial and the finite states are close to each other 
and both are ferromagnetic. As shown ini^ ior h,ho < 1 
the magnetization profile mi{t) for any finite / and t is 
non-negative. In the other domains of the quench [h^ 
and/or h is larger than 1) during relaxation mi[t) takes 
negative values, too. This type of oscillating relaxation is 
described qualitatively well by the semi-classical theory. 
The amplitude of the oscillations as well as the recurrence 
of the magnetization and the correlations are correctly 
described, but there are differences in the actual value of 
the frequencies. For quenches close to the critical point 
we expect the concept of isolated QPs to become invalid 
or at least quantitatively inaccurate due to the diverging 
correlation length either in the initial and / or final state. 

In finite system with open boundaries and for half- 
infinite systems wc find small deviations between the ex- 
act and the semi-classical results either when sites close 
to the boundaries are involved or for times t > T'poriod/2, 
when QPs refiected at the boundaries contribute to the 
magnetization or correlation reconstruction. A possible 
source for the deviations in the first case is the lack of 
translational invariance in chains with open boundaries, 
which results in spatially inhomogeneous creation prob- 
ability of QP pairs, at least close to the boundaries. The 
second kind of deviations could originate in the dynam- 
ical processes during the reflection at the open bound- 
aries, which might be more complicated than just mo- 
mentum inversion. Both effects are absent in systems 
with open boundaries, for which reason we expect our 
predictions to be accurate for all times in finite chains 
with periodic boundary conditions. 

Our semi-classical theory can be generalized to several 
directions. This theory is also valid for transverse Ising 
chains involving a sum over more ferromagnetic short- 
range interactions than only nearest neighbors, as has 
been argued for the equilibrium relaxation dynamics at 
finite temperatures by Sachdev and Youngi^-. The semi- 
classical theory should be applicable to non-integrable 
models, too, for which one has to include QP-collision 
and scattering processes. Here the quantum Boltzmann 
equation seems to be a promising approach^, as has been 
demonstrated recently for a bosonic system in2£. 



Appendix 

Here we compare our semi-classical calculation with 
the predictions by Calabrese et al^. First we recapitu- 
late the exact solution of the free fermion representation 
the Hamiltonian in Eq.(j4]) for periodic boundary condi- 
tions. In this case there are pairs of fermions with quasi- 
momenta p and — p, and in the ground state sector these 



functions: 



Vhip) 



< p < TT. Here wc define the 



leh{p) + h- cosp 



I Ship) - {h- cosp) 



and 



Up = Uha{p)uh{p) + Vha{p)vh{p) 

Vp = Uho{p)vh{p) - Vh„{p)uh{p) 



(48) 



(49) 
(50) 



in terms of which the ground state for t < (l^'o)) is 
expressed with the ground state at t > (|0)) as: 



l*o) = 



n [Up + ^Vp7^W-p\ |0) 



(51) 



Then the density of quasi-particle excitations is given by 
the noncquilibrium occupation number: 



fp - i^^olvlvpl^o) = \Vp 
This can be expressed as: 

fp = ^ [1 - cos Ap] 



(52) 



(53) 



where Ap is the difference between the Bogoliubov angles 
diagonalizing 'H(h) and 'H(/io), respectively: 



cos Ap = 



hoh — {ho + h) cosp + 1 
£/io(p)e/i(p) 



(54) 



If the difference between and h is small we obtain in 
leading order for the occupation number: 



/p = i(/i - hofsiVL^p 



(55) 



The results by Calabrese et ati^ can be formally obtained 
from our semi-classical expressions, if an effective occu- 
pation number is used. For example in Ea. (j32[) for the 
correlation length and in Eg. ([55]) for the relaxation time 
one should simply replace: 



fp 



1 



In I cos A„ 



(56) 



The semi-classical results then represent the leading term 
of the exact expressions. 

According to Calabrese et al^ there is an effective 
thermal (Gibbs) distribution or generalized Gibbs ensem- 
ble (GGE), which is obtained in integrable models by 
maximizing the entropy, while keeping the energy and 
other conservation laws fixed. This leads to an effective, 
node-dependent temperature: Tcs{p), which is given by: 



cos Ao ~ tanh 



^h{p) 
2Tcff(p) 



(57) 
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or expressed with fp we have: 
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f Ship) \ , 



(58) 



At the r.h.s. we have the Fermi distribution function 
with zero chemical potential, thus the GGE condition is 
expressed in the form, that the nonequilibrium occupa- 
tion number of the given mode is equal to its thermal 
occupation at the effective temperature: Tcsip). 
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